arXiv:1506.07855vl [cond-mat.stat-mech] 25 Jun 2015 


Stochastic models of classical particle pumps : 
Density dependence of directed current 


Debasish Chaudhuri 

Indian Institute of Technology Hyderabad, Yeddumailaram 502205, Telengana, India 
E-mail: debcOiith.ac.in 


Abstract. We present and compare different versions of a simple particle pump-model that 
describes average directed current of repulsively interacting particles in a narrow channel, due to 
time-varying local potentials. We analyze the model on discrete lattice with particle exclusion, 
using three choices of potential-dependent hopping rates that obey microscopic reversibility. 
Treating the strength of the external potential as a small parameter with respect to thermal 
energy, we present a perturbative calculation to obtain the expression for average directed 
current. This depends on driving frequency, phase, and particle density. The directed current 
vanishes as density goes to zero or close packing. For two choices of hopping rates, it reaches 
maximum at intermediate densities, while for a third choice, it shows a curious current reversal 
with increasing density. This can be interpreted in terms of a particle-hole symmetry. Stochastic 
simulations of the model show good agreement with our analytic predictions. 


1. Introduction 

In everyday life we often nse mechanical devices that ntilize oscillatory force along with a valve 
mechanism to generate unidirectional flow, e.g., a simple hand pump. In this, the valve breaks 
space-inversion symmetry, as a result breaking time-reversal symmetry, leading to unidirectional 
flow. The same basic principle of oscillatory forcing, and time-reversal symmetry breaking 
leading to directed motion, is utilized in operation of various molecular motors and ion 

pumps in biological cell, and to generate overall unidirectional flow of electrons in quantum 
pumps ig 0 0 n. However, in all these cases noise, stochastic or quantum, plays important 
role in the resultant dynamics. In molecular motors, e.g., repeated hydrolysis of ATP leads 
to a stochastic oscillatory energy input, and the intrinsic head-tail directionality of polymeric 
track on which the motors move, breaks the inversion symmetry to act like a valve [T]. Most 
theoretical studies of stochastic pumps have discussed properties of noninteracting system of 
particles, apart from a few exceptions [iniiiiiiiain!. 

In this paper, we present a stochastic particle pump model in which an external time-varying 
potential pumps energy, and a spatially varying phase factor of the oscillation breaks time- 
reversal symmetry to generate a directed current. We particularly focus on the effect of inter¬ 
particle interaction on the dynamics. This interaction may be incorporated via an exclusion 
process in a spatially discretized version of the dynamics. Two variants of this model have 
already been proposed and analyzed in some detail [HI m fT6] . Here we present a unified 
description of the model on discrete lattice. The model allows for several choices of hopping rates 
dependent on instantaneous local potential, where each choice obeys microscopic reversibility. 


We show that depending on this choice, one obtains different forms of density dependence of 
average directed current. 


2. Model 

We consider particles hopping on a ring, discretizing space into s = 1,..., N lattice sites, snch 
that the system size is L = in units of lattice spacing b. We assume that particles evolve under 
a position dependent weak oscillatory potential /?W = sin(nt + (j)s) where f3 = l/ksT with 
Boltzmann constant ks and temperature T, and (ps denote local phase factor. This potential 
drives the system out of equilibrium. If the driving frequency is slow with respect to the 
diffusion time scale l//o = a^/D where D denotes the diffusivity, the system of particles would 
come to local thermal equilibrium with the instantaneous local potential. We assume microscopic 
reversibility, i.e., the hopping rates are such that given the value of local potential at any instant 
of time the detailed balance condition, is obeyed. In this relation Us 

stands for the occupation number of s-th site and is the time-dependent hopping rate 

from s-th to s ± 1-th site. At each moment the system tries to reach equilibrium distribution 
corresponding to the instantaneous potential, but lags behind as the potential itself changes 
with time. This keeps the system out of equilibrium. Following three choices of particle hopping 
rates obey local detailed balance: 

(A) = /o exp[/3V's], a symmetric hopping rate that depends only on the on-site potential 
energy [HI US]; 

(B) = /oexp[—/3(I4±i — 14)72], depends on relative strength of the potential 
energies [T6] : 

(C) = /o exp[—/3I4±i] depends only on the potential energy at the site where the 
particle hops to. 

In this paper we present a unified derivation of the time averaged DC current obtained for 
these three cases. While models B and C are able to pump DC current even in a non-interacting 
system of particles, in model A interaction is crucial in order to achieve pumping. 

We present analytic results using a perturbation theory proposed in Ref. m, and compare 
them with numerical simulations. We present a detailed comparative analysis of the three models 
of proposed above, two of which (models A and B) were already discussed earlier [T5lfT6] . 

and the third one (model C) being the main new result of this paper. 

A hard core repulsion between particles is modeled by exclusion process, in which two particles 
can not occupy the same lattice site. With this restriction, the local density ps = (n^) and two- 
point correlation functions Cg^p = {ngUp) obey the following dynamics. 
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The last equation is for the special case of nearest neighbor correlations. Writing pg = (775) and 
the multi-point correlations as Cg^p^,,, = {ugUp ...) we can re-express the above relations as 
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Thus dynamics of each order of correlation depends on correlations of higher order, following a 
Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. Note that the evolution of local 
density may be represented as dps/dt = Js-i,s — Js,s+i where the local current 

Js—l,s — {'^S—l,SpS — l '^s,s—lps) {'^S—1,S '^S,S—l)Cs — l,S (7) 

In the time periodic steady state, local current averaged over the period r = 27r/ri is independent 
of position. Therefore the net time- and space- averaged directed current is given by 

dtJs-i,s. ( 8 ) 

s=l “'0 

For potential strength As = 0 at all lattice sites, the above model reduces to homogeneous 
symmetric exclusion process, characterized by the following local density and correlation 
functions 
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etc. m, where n is the total number of particles. As we show in the following, the BBGKY 
hierarchy separates by order, if one expands local quantities like ps, Cg^p etc. in perturbative 
expansion around As = 0. This allows one to obtain exact expressions within perturbative 
expansion. 


3. Perturbative calculation 

We consider driving at all the sites with constant potential strength and frequency 

PVg = Asin(nt + (ps) = A X Us (10) 


where 
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For small values of A we can linearize the transition rates to obtain : (A) u;s,s±i = /o(l + ^^s)) 
(B) Ws,s±i = /o{l - 5 A(us±i - Us)} (C) Ws,s±i = /o(l - Ausii) and the corresponding bond 
currents are expressed as 
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Therefore the space-time averaged directed current is 
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Note that, for non-interacting particles the correlation function Cg-i^g = 0 leads to Ja = 0. That 
means, within model-A, the pump will drive particles in an averaged unidirectional fashion only 
in the presence of interaction - free particles can not be pumped within this model |15j . However, 
for the other two models this requirement is absent. Even free particles may be pumped in a 
unidirectional manner. 

Let us write down the density and correlation functions as perturbative expansion in potential 
strength A (<C 1) 


pg = p+ ^ 

Cg,p = C(2)+ ^ A^CiJ. 


(13) 


We consider the case where the phase factor of driving potential (ps = <ps where (p = /L, 

remembering L is expressed in units of lattice parameter. Within the perturbative expansion, 
the above relations for directed current may be calculated exactly usms]. The results for the 
first two cases were derived earlier, and the third one is presented in this paper. In what follows 
we derive all three results. However, before we start deriving them, let us enlist the expressions 
for current corresponding to the three variants of the model here, 
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with Qo = p — kf) = — 

constant, qq ss p(l — p) and ko 
and Jc ~ p(l — pY (See Fig0. 


Note that in the limit of large n and L keeping p = njL 
» p2(l - p). Thus Ja ~ p^(l - p), and Jb ~ p(l - p)(l - 2p), 


3.1. Time evolution and solution 

Let us first consider model C, which constitutes the main new contribution of this paper. Using 
the perturbative expansion of Eq.(13), time evolution of the first order perturbations can be 








Figure 1. (Color online) Directed current J as a function of mean density p. The points denote 
Monte-Carlo result and the lines are plot of the functions in Eq|14[ In plotting the functions we 
replaced go by p{l — p) and fco by p^{l — p). Left panel: Shows plots of analytic functions for 
models A and C, and simulation data for model C. Right panel: Comparison of simulation results 
of model B with analytic prediction. The parameters used are system size L = 16, bare hopping 
rate due to free particle diffusion /o = 0.34, potential strength A = 0.5, frequency of oscillation 
D = 0.27r such that the time period r = 10. We choose phase difference between consecutive 
lattice points (/> = 7r/2. The data (points) were collected over lOOr after equilibration over lOOr. 
All the data for model C were averaged over 10^ initial conditions, and for model B over 10® 
initial conditions. 


written as, 


dp. 

dt 


( 1 ) 


dC: 


( 1 ) 


dC. 


dt 

( 1 ) 


dt 


fo^sP^/^ + foqo^sUs, (15) 
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where AgPs^p = gs+i,p + 9s-i,p — ‘^9s,p- Note that the above time evolution, for first order terms 
in perturbative expansion, remains the same for all the three variants of the model considered 
above. This is easy to see by comparing with Ref.s dale]. 

These linear differential equations can be solved exactly to find long time limit of time-varying 
steady state [TH] . 
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Using this in Eq.(15) we find, for all s. 
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Clearly this equation can be written in the operator form, 
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with matrix elements 
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where A satisfies the eigenvalue equation A\q) = eq\q) with €q = —2(1 — cosg') and the 
eigenfunction 'ipsig) = {s\q) = {1 /VL) ex.p{—iqs) where q = 2Trk/L = (j)k where A: = 1, 2,..., 
such that ips+N = ips- Similarly, Z\q) = (e^ — iQ/fo)\q). Thus the solution may be expressed as 
\A) =-qoZ-^A\r]). 

In the real-space representation (s|A) = —qQY^^^{s\Z~^\q){q\A\m){m\r]) = —qoYPqmi^q ~ 
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where = — 2(1 —cos(/>). The equation for two point correlation function can also be solved [16 
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where = {ko/qo){Ai^'‘ + 


3.2. Averaged directed current 

Using the above relations, one can calculate the space-time averaged directed current in the 
time-periodic steady state for all three variants of the model, through the relations shown in 
Eq.(12). For the particular case of model C, using the last relation in Eq.(12), we have 
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where in the last step we have used the fact that after integration over a period r = 27r/n, only 
the time-independent combinations of Uspi^^ terms which can be expressed as 2Re[7y*Ai^^] etc. 
remain non-zero. 

Given that ps = (—*/2)e*'^'^, and as we may write = (zg'o/2)e*'^'^a where a = 

~ *^//o]) the terms in Eq.(23) may be evaluated. We find that p^A^J^ = —q^a/A, 
and = —9oa/4. Thus the second term in the parentheses {p*Ai^'^ — = 

0. Similarly one can show that = (9o/2) sine/) (ia). Thus the terms 
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where in the last step we used the expression for In the limit of large system size 

go — fco ~ p{l — pY, and thus Jc ~ /o(l — pY- 

Similar arguments may be used to derive the results corresponding to models A and B. For 
example using relations in Eq.(12), 
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One can show that (r/*—? 7 *_]^)(Ai^^+A^^^^) = (( 70 / 2 ) sine/) (ia), and thus 2 Re (?/* — ? 7 *_i)(As^^ + 
go sin(/)Im(a). Thus 
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Again, using Eq.(12) 
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The relation (r/*-77*_]^)(A^J^+Ai = (go/2) sin 0 (ia), leading to 2 Re (p* - ? 7 *_]^)(A;_\ + 

qo sin (/> Im(a). Thus one gets 
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4. Simulation 

A detailed numerical simulation for model A was presented earlier in Ref. M- Here we perform 
Monte-Carlo simulations of the models B and C and present density dependence of directed 
current J in Fig. In the stochastic simulation, we randomly choose a lattice site s with 
uniform probability and perform a trial move with rate The trial move is accepted if 

the new site s ± 1 is empty, else it is rejected. A sweep of n trial moves for a system having 
n particles is considered as one Monte-Carlo step. We use periodic boundary condition. Note 
that in simulations we do not use the linearized versions of hopping rates, unlike in theory with 
small potential strength A. Instead we use the full non-linear forms. In all our simulations we 
keep A = 0.5, unlike the perturbation theory where we assumed A <C 1. At the time-periodic 
steady state, the current is measured on each bond and then averaged over all bonds in the 




















system, and several time-periods. We also average over many initial conditions to obtain better 
statistics. For details of the parameter values used in simulations, see figure caption of Fig. 
All the simulation data show good agreement with predictions presented in Eq.(14). 

For all the three variants of the model cnrrent vanishes as packing fraction p —)• 0 and p —)• 1 
(close pack) limits. The first vanishing is due to absence of particles to carry current, and the 
second one is due to complete jamming. If all the lattice sites are occupied, within the discrete 
lattice random sequential dynamics, particles can not move. However, the detailed density 
dependence of current J shows three very different form for the three variants of the model 
considered. While for model A, Ja ~ p^ (1 — p), model B shows a dramatic effect of current 
reversal with changing density. For model B, the density dependence is Jb p(l -p)(l - 2p), 
with a new zero in current appearing at the half filling p = 1/2. This particular model has 
a symmetry under the exchange of particles with holes together with swapping direction from 
right to left. Thus a phase factor (j) that leads to free particle motion towards right, which is 
the dominant mode at low densities, will lead to free hole motion to right at high densities. 
Therefore, particle current changes direction from near p = 0 to near p=l. Atp=l/2 the 
particle and hole cnrrents cancel each other leading to Jb = 0. We performed simulations for 
model C as well, and present the numerically obtained Jc in Fig. Onr simulation resnlts 
for model B and C agree well with theoretical predictions (see FigHH. Note that, for model 
C, theory predicts a density dependence Jc ~ p(l — p)^- In Fig. [Ipwe have also plotted the 
theoretical prediction for model A, for comparison. 


5. Outlook 

We presented a discrete pump model in which an external traveling wave potential leads to 
average directed motion of particles interacting via exclusion process. We discnssed three 
possible choices of external potential dependent hopping rates, all of which obey the microscopic 
time-reversal symmetry. We studied how a resultant directed current depends on average density 
of particles. Using a perturbative expansion for small strength of external potential with respect 
to thermal noise, we obtained analytic expressions for directed current, and compared our results 
with direct Monte-Carlo simnlations to find good agreement. While the time evolntion of first 
order perturbation in local density and correlation functions, are independent of specific choice of 
the three variants of the lattice model discnssed here, the expressions for directed current depend 
on the choice of local hopping rates. The dependence of average directed cnrrent on frequency 
and phase is the same across all the three choices of hopping rates [see Eq.(12)]. However, it 


is important to note that the detailed density dependence is very different in the three choices 
of models (hopping rates) discussed - while model A predicts Ja (1 — p), model B and C 
predict Jb ~ p(l — p)(l — 2p) and Jc ~ p(l — p)^ respectively [mean field limits of Eq.( 12)]. The 
hopping rates chosen for model-A fails to generate any directed current in absence of particle 
exclusion - density correlation turns out to be necessary. On the other hand, models B and C 
allows for driving of directed current for non-interacting particles, also. In Ref. |16) we argued 
that the model B is a natural choice, if one starts from the corresponding Langevin equation 
and discretize its dynamics. This model shows a cnrious current reversal with increasing density 
of particles. 

However, later studies by us in continuum model showed very different density 
dependence [l8], in particular, absence of current reversal predicted by the discrete exclusion 
process in model-B. In order to discnss the continuum limit of our calculation, let us now nse 
the lattice parameter b explicitly. In the continuum limit, one has to take the lattice parameter 
b/L 0, with packing fraction pb 1. Since b defines the length scale of inter-particle repulsion 
as well, in this limit, one obtains results valid for non-interacting continuum dynamics. However, 
in the real continuum system the space is continuum, but the hard core particles have finite size 
and repel each other. Thus the continuum limit of the discrete exclusion process used here, 





fails to capture the behavior of hard core particles moving in continuum space under stochastic 
thermal force, and time-oscillatory external potential. A correct discrete model would require 
the length scale of exclusion process to be defined as a new variable a = vb, such that in the 
continuum limit, 6/L —>■ 0 with z/ —>■ oo keeping a constant. Experimental realization of the 
model presented here looks possible, using colloidal particles confined in narrow channels driven 
by traveling wave potential. This would provide better insight into driven many-body dynamics, 
and could have potential applications. 
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